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ABSTRACT 

It is shown that the dissociation of molecular oxygen in the polar 
upper atmosphere caused by auroral particles is accompanied by the 
increase of ozone in the dissociating layer, contrary to that due to the 
solar UV-radiation (Maeda, Aikin, 1967). Since hydrogen-ozone reaction 
is one of the most effective sources for the atmospheric hydroxyl emis- 
sion (OH vibrational-rotation band airglow), one would expect an 
enhancement of the OH-emission, following auroral 0 2 -dissociation. 
Present calculations based on the time dependent chemical equations for 
the hydrogen-oxygen atmosphere indicates, however, that the intensity 
of OH-infrared emission does not enhance during usual quiet auroras, 
except in the lower borders of the so-called (purplish- red) Type B 
auroras. The extension of the Monte Carlo calculation of the electron- 
diffusion in the upper atmosphere (Maeda, 1965) to the case of the 
atmosphere with constant vertical magnetic field is also presented, 
indicating a drastic confining effect of the magnetic field on the hori- 
zontal spread of the impinging electron flux. 


^Maeda, K. and A. C. Aikin; Temporal Variations of Mesospheric 
Oxygen and Ozone During Auroral Events, NASA Tech. Note D-3993, 
(1967) 

(2)Maeda, K. ; Diffusion of low energy auroral electrons in the atmos- 
phere, J. Atmos. Terr. Phys. 27^ 259-275, (1965) 
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1. Introduction 


The rate of dissociation of molecular oxygen by auroral electrons 
and consequent variations in the temporal and the vertical distributions 
of the oxygen allotropes were calculated with respect to the isothermal, 
pure -oxygen atmosphere (Maeda and Aikin, 1967 a, b*). 

In this paper, similar calculations are made with respect to the 
non-isothermal, hydrogen-oxygen atmosphere to see the subsequent 
auroral effects on the hydroxyl emission. 

The previously reported Monte -Carlo calculation of electron dif- 
fusion in the upper atmosphere (Maeda, 1965 a, b) is also extended to 
the atmosphere with a constant vertical magnetic field. Although this 
new result of the diffusion of auroral electrons in the polar atmosphere 
does not alter the previous conclusions (MA-I and II), the effect of the 
geomagnetic field on the horizontal spreads of impinging electron beams 
is remarkable. Since these are not essential for the present problem, 
the details of the computations on the electron diffusion in the atmosphere 
with the geomagnetic field will be presented elsewhere (Berger, Maeda 
and Seltzer, 1968). 

2. Auroral Dissociation of the Atmospheric Oxygen 

2.1 Dissociation of Oxygen Molecules by Electron Impacts. 

Due to the great improvements in the energy-resolution and con- 
trol of mono-energetic electron beam in the laboratory experiments, 
the cross-sections of molecular dissociations as well as excitations 
and of ionizations have recently been measured very accurately, almost 
comparable to those of optical experiments. 

The oxygen molecule has the large electron-impact cross-section 
for dissociation with a maximum around 8.4 eV, corresponding to the 
photodissociation due to the Schumann- Rung e UV radiations (1400 
- 2000 A). 

These cross-sections are shown in Figure 1 as a function of the 
energy of colliding electrons. In the figure, the meanings of each 


*These will be referred to as MA-I and MA-II, respectively. 
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cross-section are as follows; (i) cr“> the dissociative attachment of an 
electron i.e. 


0 2 + e '► O 4 O 


( 2 . 1 ) 


This has a peak (1.6 x 10 18 cm 2 ) at 6.5 eV (Rapp et al., 1964 a, b). 
(ii) a., the total ionization, including dissociative ionization such as 


0 2 4 e — > Q 4 0 + 4 2e 

O 4 0 + 4 e (2.2) 

O' 4 0 + 4 3e etc. 


With increasing energy of electron, the multiple ionization (0 ++ , 0 + + + ) 
appears in the secondary products (Frost and McDowell, 1958). In this 
figure, o', is taken from the result given by Tates and Smith (1932). 

(iii) o SL ; molecular dissociation into two neutral atoms, i.e. 


0 2 4 e 4 0* 4 e 


(2.3) 


where O and O* are mostly 3 P and ! D-states, respectively. It should be 
noted that this cross section is the sum of several cross.-sections cor- 
responding to the various transitions between the initial states of oxygen 
molecules (not only ground state X 3 2 g but also a and b ) and 

the final dissociative states (such as B 3 ^ u ►~0( 3 P) 4 0( 1 D) # 

1 A y ►0( 1 D) + O(’D) and ^ *-0( 1 D) + 0( 1 S) etc.). 

According to Lassettre et al. (1964), and Silverman et al. (1964), 
the maximum cross-section for the electron impact 0 2 -dis sociation 
is of the order of 4 x 1 0 16 cm 2 above 8.4 eV. It should also be noted 
that above approximately 18.7 eV dissociation is accompanied by 
ionization as shown in (2.2). 

the cross-section corresponding to (2.3) calculated by 
Bauer and Bartky (1965) based on the classical theory of Gryzinski 
(1959), which gives a lower bound on the electron impact dissociation 
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cross section, since only one transition has been taken into account. 
cr M is the similar cross section used by Maeda (1962, 1963) which was 
derived by combining the measurements of Glockler and Wilson (1932) 
below 20 eV and theoretical extrapolation above 100 eV (Massey and 
Burhop, 1956). Figure 1 indicates that the rate of 0 2 -dis sociation by 
electron impacts based on cr M was nearly two order of magnitude under- 
estimated as compared to the latest laboratory data. 

Finally, the recent semi-empirical expression of electron impact 
cross-sections for molecular oxygen given by Green and Dutta (1967) 
and Watson et al. (1967) agrees with cr. shown in Figure 1. Their 
cross section for excitations into the Schumann- Rung e band and con- 
tinuum seems, however, to be too small, if not neglected, even compared 
to theoretical results of Bauer and Bartkey (1965). 

2.2 Diffusion of Auroral Electrons in the Polar Atmosphere. 

Due to the fluctuations in energy loss (straggling) and multiple 
scattering, the energy distribution of penetrating electrons in the at- 
mosphere becomes broader with increasing depth of penetration even 
if the initial primary electrons are monoenergetic. 

This is shown in Figure 2 where i(E 0 , E, ^); the number of elec- 
trons with energies between E and E + dE at the depth correspond- 
ing to the incident energy E 0 , is plotted against E with the parameter 
Non-dimensional penetration depth g is defined by the ratio of the 
actual penetration depth in the atmosphere (in g/cm 2 ) to an empirical 
range r 0 (in g/cm 2 ) which is given by 


r o = cE o 7/4 (2.4) 

where E 0 is the initial energy in kev and c = 4.57 x ICf 6 (Grun, 1957). 
It should be noted that these calculations are possible only by means 
of the Monte-Carlo method (Berger, 1963). 

For the practical applications, the results of Monte-Carlo calcu- 
lations are approximated by the empirical formula (MA-I, II). These 
are also shown in Figure 2 by dashed lines. As already mentioned, 
in these calculations of electron diffusion in air, the energy losses by 
atomic collisions, multiple coulomb scattering as well as straggling 
(statistical fluctuations in energy loss) are simultaneously taken into 
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account, but the effect of the earth's magnetic field is, however, com- 
pletely neglected. It was shown, however, that as long as the energy 
dissipations in each horizontal layers are concerned, the vertical mag- 
netic field parallel to the incident direction has no effect on the results 
(MA-I, II). The effect of the vertical magnetic field is simply to re- 
duce the horizontal spread of electron beam by curving the orbits of 
eachfree electrons between their collisions in air. This effect is, 
however, very remarkable as shown in Figure 3 where the mean hori- 
zontal radius, R 0 , of vertically incident electron beam with initial 
energy E 0 = 50 kev is plotted against altitude z (in km), for three 
cases of vertical magnetic field H i.e. H = 0, 0.3 and 0.6 gauss. The 
definition of mean horizontal radius R 0 and some other related quantities 
are discussed in Appendix A. 

From this figure, one can see that the horizontal spread of auroral 
electron beams are confined by the vertical earth's magnetic field of 
0.6 gauss within the range of less than 100 m, which should extend 
upto nearly 10 km, if there is no geomagnetic field. 

In Figure 4, the gyroradius R c (in meter) of 50 kev electron in the 
0.6 gauss field and the gas -kinematic mean free path R f (in meter) for 
50 kev electron are shown for comparison, indicating the importance 
of magnetic control of horizontal spread even below the altitude of 
100 km, where collisions dominate i.e. the collision frequency between 
electron and air molecules and ions, v e# exceeds the gyrof requency of 
electron in 0.6 gauss, ^ c (Chapman, 1956). 

Finally, the horizontal spread of auroral electrons and their verti- 
cal distributions are shown in Figure 5 for incident energies, E 0 = 400/n 
kev, where n = 1, 2 • • • 8. The results for E 0 < 4 kev are, however, 
not conclusive because of inaccuracies in energy loss near the terminat- 
ing point of electron path i.e. the relative accuracy decreases near the 
cut-off energy computation (Berger, Maeda and Seltzer, 1968). 


2.3 Rate Coefficient for Auroral 0 2 -Dissociation 

The rate of dissociation of atmospheric oxygen molecules by the 
monoenergetic primary electron, E 0 , at the altitude z can be given by 


Jde ( Z ’ E 0) 


f E ° 

~f E 

N(z) 

CT d O) j e O’ E > z ) dw 

'0 

L_ c _J 


i(E 0 , E, z) dEdz (2.5) 
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where 


N(z) is the number density of molecular oxygen at the altitude z, 
(cm ~ 3 ). 

o- d (w) is the differential cross-section of O a -dissociation by elec- 
trons with energy w (cm 2 , w is kev). 

j e (w, E, z) is the number of secondary and tertiary electrons of 
energy between w and w + dw, produced by electrons of energy E per 
unit thickness of air at the altitude z(cm _1 , z in km). 

i(E 0 , E, z) is the number of electrons with energies between E 
and E + dE at the altitude z, corresponding to the incident energy E 0 . 

E c is the threshold energy for 0 2 -dissociation and assumed to be 
5 eV. 

As mentioned in Appendix A, transformation from i(E Q , E, g) to 
i(E 0 , E, z) is made with respect to the CIRA (1964) atmosphere. The 
actual calculations of J de (E 0 , z) are performed numerically, using 
o- gL shown in Figure 1 as cr d (w) with several approximations (MA-I 
and II). The results are shown in Figure 6. 

According to the direct measurements by rockets (Mcllwain, I960; 
Evans, 1966) and by satellites (O'Brien, 1964), the energy spectrum as 
well as intensities of auroral electrons fluctuate widely within the 
order of hundred meters and within afraction of a second. If one take s the 
spatial and temporal average over the range of several hundred meters 
and several minutes, the incident auroral electron spectra are ex- 
pressed with a good approximation by 


i(E 0 ) = 

where E is a constant and around 10 kev or less in most cases of 

a 

bright aurora. E q is reported occasionally as high as 100 kev during 
the events of auroral absorption (Brown, 1964, 1966; Bailey et al., 1966). 
The rate of 0 2 -dissociation by the auroral electrons with these energy 
spectra are calculated as follows: using the previous results for mono- 
energetic electrons; 


i n exp 


(T) 


( 2 . 6 ) 


5 


f 

J n 


” E n /E 

T (z, E 0 ) c 0/ a dE r 


J D (z-E a ) = 


r 

J n 


exp 


dE r 




i v * 1 

En Jde (Z ’ E °‘ )eXp 

i = 1 



where J de (z, E Q A ) r s are given by (2.5) and shown in Figure 6. 


(2.7) 


Considering two rather extreme cases of soft and hard spectrum 
of auroral electrons, J D (z, E a )'s are calculated for E a = 5 kev (soft 
spectrum) and 100 kev (hard spectrum) as shown in Figure 7. 


3. Variations of Oxygens and Hydroxyl During Auroral Events 


3.1 Equations and Constants for Chemical Reations. 


All reactions in hydrogen-oxygen atmosphere are discussed by 
Hesstvedt (1964). The principal reactions considered in the present 
calculation are shown in Table 1. In this table, J 1? J 2 * * * J 4 and k lf 
k 2 ° * k 20 are the rate coefficients for dissociations and each chem- 

ical reactions. The numerical values of these coefficients and their 
temperature dependences are shown in Table 2. Some of these coeffi- 
cients are not known accurately as others such as , k 2 and J 2 . 

3.2 The First Approximations (Constant Nighttime Atomic Hydrogen). 

Since the main purpose of the present paper is to see the effect of 
auroral electrons on the hydroxyl night sky emission, the following 
three equations are used in the preliminary analysis. This approxima- 
tion was used by several investigators for the nighttime (after sun-set) 
enhancement of hydroxyl emission (Wallace, 1962; Ballif and Venka- 
teswaran, 1 963) 

= 2J 2 t0 2 J + J 3 [0 3 ] - 2k j [M] [O] 2 - k 2 [0] [0 2 ] [M] 

- k 3 [0] [0 3 ] " k 6 [01 [OH] - k 8 [0] [H0 2 1 (3.1) 

dLo 3 ] 

r - k 2 [0] [0 2 ] IM] - k 3 [0] [0 3 ] - k 4 [H] [0 3 ] - J 3 [0 3 ] (3.2) 
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d [H] 


dT = [H 2 0] + kg [OH] [0] - k 4 [H] [0 3 ] - k 9 [H] [0 2 ] [M] (3.3) 


[0] + 2 [0„] + 3[0,] = f(z) 


(3.4) 


where J = + J g and f(z) is a function of altitude. During the night 

time, photo-dissociation is zero (i.e. = 0) and variations of hydro- 

gen atom can be neglected (Wallace, 1962) then from Equation (3.3), 
we get 


[OH] = [H] (k 4 [0 3 ] +k 9 [0 2 ] [M])/k 6 [0] (3.5) 


Since it is shown that the soft spectrum auroral electrons cause no 
significant effects on the oxygen allotropes concentrations in the upper 
atmosphere (MA-I, II), only the hard spectrum case i.e. E a = 100 kev 
in Equation (2.6) is considered in this paper. 

The ineffectiveness of soft auroral electrons on the atmospheric 
oxygen allotrope variation is due to the very large time constant (i.e. the 
life time) of atomic oxygen above 90 km level, where most of the oxygen 
molecules are dissociated by the solar UV-radiation. The decay time 
r x (in sec.) of the mesospheric oxygen atom is shown in Figure 8, 
together with that of ozone r 3 (in sec.). It should be noted that the 
life time of atmospheric ozone decreases slowly with altitude, while 
that of atomic oxygen increases drastically above 90 km. The similar 
curves given by Hunt (1965) are also shown in Figure 8 for a comparison. 

Temporal variations of the vertical distribution of the oxygen atom 
and of ozone caused by the hard spectrum auroral electrons are shown 
in Figure 9, where the flux of the primary auroral electron is assumed 
to be 10 10 electrons per cm 2 sec. The corresponding time variations 
of hydroxyl concentration calculated under the constant nighttime con- 
centration of atomic hydrogen are shown by dashed lines in Figure 11. 

As shown by Hunt (1966), the daytime concentration of atomic hydrogen 
increases below around 85 km, which is indicated by the dashed line in 
Figure 9. 

The result which is based on the time independent nighttime dis- 
tribution of atomic hydrogen shows that the hydroxyl concentration does 
not inhance but rather diminishes during hard spectrum auroral events. 
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The OH-emission corresponding to the initial condition {t = 0 min) 
indicated by the dashed line in Figure 11, is of the order of 2000 R, 
provided that the F-factor* is assumed to be 4 (Chamberlain and Smith, 
1959). 


3.3 The Second Approximation (Time Dependent Nighttime Hydrogen). 

Considering all chemical reactions listed in Table 1, we get eight 
simultaneous time dependent differential equations for the chemical 
variations at each levels in the atmosphere as shown in Appendix B. 

Two differential equations are, however, replaced by the additional two 
conditions which indicate the conservations of oxygens and hydrogens 
at each level, respectively. These are also shown in Appendix B. Actual 
calculations to see the auroral effects are made after finishing the 
marching integration for a period of five days. In other words, the 
deviations due to the arbitrary initial conditions (i. e. H0 2 , HO, O and 
0 3 an be assumed to be zero at the beginning, for example) are elimi- 
nated by computing the diurnal variations of these compounds, including 
the photodissociations by the solar UV-radiations. 

The variations of the concentrations of atomic oxygen [o] , ozone 
[Oj and atomic hydrogen [h] due to the auroral electrons with flux 10 
electrons per cm 2 per second are shown in Figure 10, and the corre- 
sponding variation of OH is shown in Figure 1 1 by full lines. 

4. Conclusions 

The vibration- rotation bands of hydroxyl (OH) discovered by Meinel 
(1950) is one of the brightest airglow in night sky. The enhancement 
of this emission after the sunset, which sometimes exceeds 5000 kR, 
is well -explained by hydrogen-ozone reactions (2.4) in Table 1 pro- 
posed by Bates-Nicolet (1950). As shown by Wallace (1962), so far the 
following observational features have been known: 

(1) The total intensity of emission at its maximum which occurs 
around hours after the sunset is of the order of 5000 kR or 
more. 

(2) The height of emission is between 75 km and 90 km with a 
maximum around 85 km. 

* Average number of photons emitted per excited OH -molecule (cf. Hunt, 
1966). 
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(3) The emission ranges from hardly visible 5000 A region to 
far infrared 4.5 /z with maximum around Z.3 ~ 3.6 (jl . 

(4) The emission band intensity drops sharply beyond v = 9, indi- 
cating that the excitation of atmosphere hydroxyl is up to 3.3 
eV. 

(5) Generally, the emission is stronger in winter than in summer, 
although day-to-day fluctuation is occasionally very large. 

(6) Relations of OH-emission intensity and its fluctuation to the 
geomagnetic activity and solar activities are not well-estab- 
lished. 

Considering the above feature, especially the nighttime enhance- 
ment due to hydrogen ozone reaction, one would expect the enhancement 
of OH-emission during hard spectrum auroras. 

The present calculation shows, however, that hydroxyl emission 
does not enhance but rather decreases, if the nighttime concentration 
of hydrogen atom is assumed to be constant as shown by dashed lines 
in Figure 11. This indicates a fundamental difference between the photo- 
dissociations of the atmospheric oxygens due to the solar UV-radiation 
and the auroral dissociations. This difference can be attributed to the 
steep decrease of the solar UV-intensity in the shorter wavelength. 
Namely, the intensity of O 2 -dissociating o UV-radiation, i. e. , Schumann - 
Runge band and continuum (1400 ^2000 A) is nearly 10” 4 times of O s - 
dissociating radiation, (i.e. , Hartley band, 2200 ~ 3000 A and some 
visible radiations in Chappins band 4500 ~ 7500 A)* * 

In the case of the photo-dissociation by the solar UV-radiation, the 
increase of atomic oxygen is therefore always accompanied by the de- 
crease of ozone, since the photo-dissociation of ozone (the last term 
in the right side of Eq. (3.2)) dominates over the production by three 
body collision (the first term in the right side of Eq. (3.2)). 


*The solar spectral irradiances are roughly 10 2 (ergs /cm 2 sec m) at 
1500A, 3.10 5 (ergs/cm 2 sec M) at 3000A and 2.10 6 (ergs/cm 2 sec M) 

o 

at 5000A, respectively. 


9 


On the other hand, in the case of auroral dissociation, even the 
cross-section of ozone is assumed to be more than one hundred times 
that of molecular oxygen; 0 3 -dissociation does not dominate over the 
production by three body collisions below around 80 km. Consequently, 
an increase of atomic oxygen caused by O 2 -dissociation is associated 
with increase of ozone. 

As can be seen from Equation (3.5), the time variation of hydroxyl 
concentration is essentially proportional to the ratio of atmospheric 
ozone concentration to oxygen atoms concentration if the variation of 
hydrogen atom is ignored. Consequently, the concentration of hydroxyl 
does not increase during auroral events, but rather decreases. 

On the other hand, all possible chemical reactions between oxygens 
and hydrogens are considered, the concentration of atomic hydrogen 
below 80 km level increases during the auroral events as shown in 
Figure 10, even though the direct dissociations of water vapor by auroral 
electrons are neglected. Correspondingly, the concentration of OH in- 
creases as shown by full lines in Figure 11. The altitude of the maximum 
OH-concentration is, however, lower and its vertical distribution is 
broader than the observed concentrations of OH-emission (Packer, 1961). 
These results indicate that the auroral enhancement of OH-emission 
should decay faster than the usual (quiet time) nighttime OH-emissions. 
Although it is not shown in the present figures, the enhancement of OH- 
emission based on the second order approximation continues several 
minutes after the cessation of electron precipitations and decays quickly 
in the lower atmosphere. The maximum OH-emission expected from 
the present results (the second order approximation) is only around 
400 R (with F-factor 4). If the direct auroral H 2 0-dissociation is 
taken into account, the maximum enhancement of OH-emission should 
exceed far beyond this value. In the case of strong auroral glow with 
very large flux of high energy protons, another additional enhancement 
of OH-emission can be expected together with well-known hydrogen 
emissions. It should be noted that if the flux of auroral electrons is 
less than 10 5 per cm 2 sec, no auroral effect can be observed. This can 
be explained by the domination of loss term over production term in the 
time dependent chemical equations for atomic oxygen and for ozone. 

In the present calculation, the vertical temperature distribution 
which affects on the rate coefficients as shown in Table 2, is taken 
from CIRA (1964) atmospheric model. Although it is not expressed 
explicitly in the present calculations, the photodissociation of water 
vapor and methane are regarded as the sources of hydrogen atoms in 
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the upper mesosphere (Nicolet, 1965). It should be noted that in the 
altitude range where the variations of all oxygen allotropes and hydrogen- 
oxygen compounds take place, i. e. above 50 km up to 90 km, the dy- 
namical effects of atmospheric motion especially diffusions due to 
turbulence cannot be neglected. In this respect, the present calculations 
need further improvement, including the auroral dissociations of water 
vapor and contributions from auroral protons. 

As a summary, the present calculations (the second order approx- 
imation) indicate that the hydroxyl infrared emission will be enhanced 
during strong Type B auroras (the so-called purplish-red lower border), 
due to the auroral dissociation of atmospheric oxygen molecules and 
subsequent increase of ozone and atomic hydrogen reaction in the 
upper mesosphere. 

Table 1 

Principal Reactions in Hydrogen Oxygen Atmosphere 






Table 1 (Continued) 


(2) Chemical Reactions by Collision 

k. 

0 + 0 + M ^Oj + M 


O + 0 2 + M- 


-0 3 + M 


o + o 


-o 2 + o 2 


H + 0 3 - 


•OH + O 2 


■ho 2 + O 


OH + O: 


HO 2 + O- 


k, 

k: 


:H + 0 2 


-OH + 0 2 


H + 0 2 + M- 


-H0 2 + M 


H + O + M ►■OH + M 


H + H0 2 - ^ ■■■-■-*>H 9 0 + O 
H 2 + 0 2 


OH + OH 



H + H + M- 


-^-►H 2 + M 


H + OH , 1 5 te H 0 + O 


OH + O, 


I 

k 


H 2 + OH 2Z_*H 2 0 + H 


il*.ho 2 + o 


OH + OH ^H, O + O 


OH + H0 2 ►H 2 O + 0 2 


( 2 . 1 ) 

( 2 . 2 ) 

(2.3) 

(2.4) 

(2.5) 

( 2 . 6 ) 

(2.7) 

( 2 . 8 ) 

(2.9) 

( 2 . 10 ) 

( 2 . 11 ) 

( 2 . 12 ) 

(2.13) 

(2.14) 

(2.15) 

(2.16) 

(2.17) 

(2.18) 

(2.19) 

( 2 . 20 ) 



Table 2 

Numerical Values of Chemical Rate Coefficients 
Corresponding to the Reactions Shown in Table 1 

kj = 5 x 10 " 34 /T (cm 6 sec 3 ) (Nicolet-Bates, 1950) 

k 2 = 8 x 10 " 35 exp(-^r-) (cm 6 sec” 1 ) (Benson-Axworthy, 1965) 

C Z S 5 0 \ 

j (cm 3 sec’ 1 ) (Benson-Axworthy, 

1965) 

/ 1500 \ , 

k 4 = 1.5 x 10 11 VT exp — y~ ) (c ™ 3 sec ) (Wallace, 1962) 

k 5 = 1.5 x 10 ” 11 fr exp(-^^) (cm 3 sec’ 1 ) (Wallace, 1962) 

k 6 = 6.7 X 10 ” 12 YT exp (- -tf-) ( cm3 sec’ 1 ) (Hesstvedt, 1964) 

k 7 = 3.6 x 10 " 9 YT exp(- 12, °^ — ) (cm 3 sec" 1 ) (Hesstvedt, 1964) 

k a = 6 x 10 ’ 13 Vt (cm 3 sec’ 1 ) (Hesstvedt, 1964) 

k 9 = 6 x 10’ 37 YT (cm 6 sec ) (Hesstvedt, 1964) 

k 10 = 6 x 10’ 34 YT (cm 6 sec" 1 ) (Wallace, 1962 ) 

k n = 1.5 x 10" 11 YT exp (- (cm 3 sec' 1 ) (Hesstvedt, 1964) 

k 12 = 2 x 10 " 13 YT (cm 3 sec" 1 ) (Hesstvedt, 1964) 

k 13 = 6 x 10 13 yT (cm 3 sec X ) (Hesstvedt, 1964) 

k 14 = 1.6 x 10 ~ 33 yT (cm 6 sec" 1 ) (Hesstvedt, 1964) 

-i! ( 2000\ 3 

k = k 5 = 1.5 xio /T exp j (cm sec ) (W allace, 

1962) 
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Table 2 (Continued) 

k 16 = 4.5 xICT 11 Yt exp ) ( cm3 sec l ) (Wallace, 1962) 

k 17 = 7 x 10 12 /T exp ^ T~') (cm 3 sec *) (Hesstvedt, 1964) 

k 18 = 1.5 x 10 11 Yt exp ( (cm 3 sec *) (Hesstvedt, 

1964) 

k 19 = 1.6 x 10 13 /T (cm 3 sec *) (Hesstvedt, 1964) 

k 2 o = k 13 = 1.6 x 10 33 /T" (cm 3 sec *) (Hesstvedt, 1964) 
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Figure 1. Total cross-sections of the 0 2 -dissociation by electron im- 
pacts. cr~ , a SL and cr. are obtained by the laboratory experiments of 
Rapp et al (1964), Lassettre et al (1964) and Tates-Smith (1932), 
respectively. cr BB and cr M are the similar cross-section as a sL* 
calculated theoretically by Bauer and Bartky (1965), and the one 
derived by Maeda (1962, 1963), respectively. 
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Figure 2. The changes of 
differential energy spectrum 
of the monoenergetically 
incident electrons, i (E Q , E, £) 
with penetrating depth g . 

The full lines and dashed 
lines correspond to the di- 
rect results of Monte Carlo 
calculation and to those given 
by an empirical formula 
(Maeda and Aikin, 1967). 
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Figure 3. Horizontal spreads of vertically incident 50 kev electrons 

in the non-geomagnetic field atmosphere and in the atmosphere with 
vertical field 0.3 and 0.6 gauss, calculated by the Monte Carlo method. 
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Figure 4. The gyroradius R c (in meter) of 50 kev electrons in the ver- 
tical magnetic field of 0.6 gauss and the gas kinematic means free 
path of 50 kev electrons R f (in meter) in the CIRA (1964) atmos- 
phere, in comparison with the present result of the horizontal spread of 
50 kev electrons in the CIRA (1964) atmosphere with 0.6 gauss field. 
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Figure 5. The horizontal spreads of electron beam in the CIRA (1964) 
atmosphere with vertical 0.6 gauss magnetic field for the incident 
energies, E = 400/n kev, where n = 1.2 ... 8 
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Figure 8. The 
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by the hard spectrum auroral electrons of 10 10 cm“ 2 sec” 1 flux in the atmosphere with constant 
atomic hydrogen distribution. The altitude distributions of atomic hydrogen [H] are shown by 
thin lines, and dashed line below around 85 km corresponds to the daytime maximum (Hunt, 1966). 
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Figure 10. The altitude dependence of temporal variations of oxygen atom [O] , ozone [0 3 ] and 
hydrogen ktom [H] caused by the hard spectrum auroral electrons of 10 10 cm 2 sec" 1 flux. 









APPENDIX A 


The mean horizontal radius R Q is defined by 




r 

r\ 


RF(E 0 . 


z, R.) dR 



z, R) dR 


(A. 1 ) 


where F(E Q , z, R) is the horizontal distribution function of energy dissi- 
pation per unit depth in a ring between radius R and R + dR, at the alti- 
tude z, corresponding to the initial (incident) energy of electrons E Q . 

In other words, if the energy deposition per unit depth at the altitude z 
is D(E 0 , z), then 


D(E 0 , 


z) = 



z, R) dR 


(in eV/g cm -2 ) (A. 2) 


the energy dissipated per unit volume, d(E Q , z) (in eV per cm 3 ) is 


d(E 0 , z) 


F(E 0 , z, R) • p(z) 

2 77 R 


(A. 3) 


— 3 

where p (z) is the atmospheric density at the altitude z (g cm" , zinkm). 
The original Monte-Carlo calculations are made, using the atmospheric 
depth x (in g/cm 2 ). Transformations of all quantities as a function of 
altitude z are made with respect to the CIRA (1964)-atmosphere. 
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APPENDIX B 


Time dependent chemical equations used in the calculations 
second order approximation are as follow; 


d[Q] 

dt 


Aj - Bj [0] - C x [0] 2 


where 

Aj = 2J 2 [0 2 ] + J 3 [0 3 ] + J 4 [H0 2 ] 

+ k 5 [H] [0 3 ] + k u [H] [HO a ] + k is [H] [OH] 

+ k lg [OH] [0 3 ] + k, g [OH] 2 

Bj = k 2 [0 2 ] [M] + k 3 [0 3 ] + k 6 [OH] 

+ k g [H0 2 ] + k 10 [H] [M] 

Cj = k t M 

d[0 2 ] 

dt _ ^ 1 ~ 

where 

A 2 = kj [o] 2 [M.l + 2 • k A fOj [0 3 j + k 8 [Hi [0 3 3 

+ k 6 [OH I [0] + k 8 [H0 3 ] [0] * k 12 [H] [H0 2 J 
+ k 20 [OH] [H0 2 ] 

B 2 -- k 2 [O] [Mj + k 7 [H] - k,, [HJ -.Ml 


of the 
(B.l) 

(B. 1* ) 

(b.l") 

(B.l"') 

(B.2 ) 

(B.2-) 

(B.2") 
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(B.3) 


d [O 3 ] 

dt A 3 ~ B 3 ^® 3 ^ 

where 

a 3 = k 2 [ 0 ] [ 0 2 ] CM] 

b 3 = J 3 + k 3 [0] + (k 4 + k 5 ) [H] + k 18 [OH] 

d[H] . r . r .. 

"dF = A 4 “ B 4 M - C 4 [H] 2 


(B.3') 

(B.3") 

(B.4) 


where 

A 4 = Jj [H 2 0] + k 6 [OH] [0] + k 17 [Hj] [OH] (B.4' ) 

B 4 = (k 4 +k 5 ) [0 3 ] +k ? fo 2 ] +k g [0 2 ] [Mfl +k 1Q fo] [M] + (k u +k 12 +k 13 ) lH0 2 ] +k, s [OH] 

(B.4") 

C 4 = k 14 [M] (B.4"') 

d[H 2 ] 

dt -A 5 _B S^ B 2^ (B. 5) 


whe re 


A 

B 


5 


5 


k 12 [H] [H0 2 ] 
k 17 [OH] 
d[ 0 «] 

dt - A 6 


+ k 14 [H] 2 [M] + k 15 [H] [OH] 


- B 6 [OH] - C 6 [OH] 2 


(B. 5' ) 
(B.5") 

(B. 6 ) 
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where 


A e = Ji [H 2 0] + J 4 [H0 2 ] + k 4 M [0 3 ] + k g [H0 ? ] [0] 

+ k 10 [H] [0] [M] + 2k j 3 [H] [H0 2 ] + k 16 [H 2 ] [0] (B.6') 

B g = k 6 [0] + k ls [H] + k 17 [H 2 ] + k lg [O 3 ] + k 20 [HD 2 ] (B. 6 ") 

C 6 = k 19 (B. 6 »«) 

d[H 2 0] 

— dT“ = A 7 - B 7 C H 2 0] (B.7) 

where 

A 7 = k M [H] [H0 2 ] + k 17 CH 2 ] [OH] + k 19 [OH ] 2 + k 2Q [OH] [HO,,] (B.7') 

B 7 = J t (B.7" ) 

d [HO,] 

dt - = *8 - B 8 lH ° 2 ] < B - 8) 

where 

A 8 = k s [H] [0 3 ] + k 9 [H] [0 2 ] [M] + k lg [OH] [0 3 ] (B. 8 ' ) 

Bg = k 8 [0] + (k n +k 12 +k 13 ) [H] + k 20 [OH] (B. 8 ") 

The numerical values of rate coefficients are shown in Table 2. 
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In the actual calculations, two equations (B.2) and (B.8) are re- 
placed by the following two conditions; 


tO] + 2t0 2 ] + 3 [O3] + 2 [H0 2 ] + [OH] + [H 2 0] = a[M] (B.9) 

[H] + 2tH 2 ] + 2[H 2 0] + [H0 2 ] + tOH] = /3[M] (B.10) 

where a and /3 are assumed to be 0.4 and 5. 10 -s , indicating a relative 
abundance of oxygen and hydrogen at each level in the atmosphere, 
respectively (cf. Hesstvedt, 1964). 
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